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ABSTRACT 

We determine the importance of redshift-dependent systematic effects in the determination of stel- 
lar masses from broad band spectral energy distributions (SEDs), using high quality kinematic and 
photometric data of early- type galaxies at z ~ 1 and z ^ 0. We find that photometric masses of z ~ 1 
galaxies can be systematically different, by up to a factor of 2, from photometric masses of z ~ 
galaxies with the same dynamical mass. The magnitude of this bias depends on the choice of stellar 
population synthesis model and the rest-frame wavelength range used in the fits. The best result, 
i.e., without significant bias, is obtained when rest-frame optical SEDs are fitted with models from 
Bruzual & Chariot (2003). When the SEDs are extended to the rest-frame near-IR, a bias is intro- 
duced: photometric masses of the z ~ 1 galaxies increase by a factor of 2 relative to the photometric 
masses of the z ~ galaxies. When we use the Maraston (2005) models, the photometric masses of 
the z ~ 1 galaxies are low relative to the photometric masses of the z ~ galaxies by a factor of 
~ 1.8. This offset occurs both for fits based on rest-frame optical SEDs, and fits based on rest-frame 
optical+near-IR SEDs. The results indicate that model uncertainties produce uncertainties as high 
as a factor of 2.5 in mass estimates from rest-frame near-IR photometry, independent of uncertainties 
due to unknown star formation histories. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Galaxy masses are an essential link between theories 
for galaxy formation and observations of the galaxy pop- 
ulation and the evolution thereof. In the local uni- 
verse masses can be measured accurately by modeling 
the luminosity distribution and the dynamical struc- 
ture of galaxies (see, e.g.. iCappellari et al.ll2005[) . Scal- 
ing relations s uch as the fundamental plane for early- 
type galaxies ijDiorgovski fc Davisl Il987t iDressler et al.l 
119871) and the Tully- Fisher relation for spiral galaxies 
llTullv fc Fish^ll97l can be used to measure the evolu- 
tion of the mas s-to-light rati o (M/L) from high redshift 
to the present l)Franxl ll993). This technique has been 
applied successfully and has provided constraints on the 
formation ep och of massive early-type gala xies out to 
z ~ 1.3 (e.g.,E Dokkum fc Staufordll200l. 

The number of galaxies with dynamically measured 
masses at intermediate redshifts (z ~ 1) is small, be- 
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cause obtaining those is observationally expensive. Fur- 
thermore, thos e samples are severely ha mpered by selec- 
tion effects (see lvan der Wel et al.ll2005j) . Therefore, it is 
not yet possible to obtain directly the galaxy mass den- 
sity at high redshift. Also, even though the redshift at 
which dynamical masses can be measured is steadily in- 
creasing, the most active era of galaxy formation, z > 2, 
is not accessible in that respect with the current gener- 
ation of instruments. For all of these reasons, one has 
to rely on less accurate mass estimates to construct a 
picture of the high-z galaxy population lShai)kiy et_al| 
| 2003HPapovich et al.l200llForster Schreiber et al.l2004t 
ILahhe et al .1120051 IShanlev et al .1120051) and the evolu- 
tion of the mass density with redshift ((Bell et al .1120041: 



iDrorv et al~ll2004al iFaber et al.ll2005|) . In these stud- 
ies, broadband photometry is compared with predic- 
tions from stellar popul ation models (most commonly, 
iBruzual fc Charlotl2003|) in order to constrain the physi- 
cal properties of high z-galaxies, and thereby their stellar 
masses. 

However, the uncertainties in mass estimates obtained 
by fitting the spectral energy distributions (SEDs) of 
galaxies are large. It is well-known that parameters 
such as age, dust content and metallicity are degenerate, 
leaving M/L uncertain. The lack of knowledge of the 
stellar initial mass function (IMF) leads to a (probably 
large) systematic uncertainty. Comparisons of photomet- 
ric mass estimates with dynamical mass measurements 
are essential to establ ish their robustness and accuracy. 
iBell fc de .Tonal 1(200 ll) have shown that the optical col- 
ors of spiral galaxies cor r elate well with their M/L. At 
z ~ 1 Ivan der Wel et aTl (|2005^l have demonstrated that 
optical colors correlate well with dynamically determined 
M/L for early-type galaxies. These results suggest that 
full SEDs should provide good constraints on masses and 
M/L. This has been tested directly for local galax- 
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Note. — Photometry of the high-rcdshift field galaxy sample. IDs are the same as in van der Wcl ct al. (2005b). All magnitudes within a 
5"-diameter aperture of PSF-matchcd images. The two lower-limits that occur in the table are 3er-limits. Object CDFS-29 falls outside the 6435 
ACS mosaic. The photometry in the columns without listed errors have errors on all individual objects of 0.05 mag or less. The typical errors on 
the 2775 , Z850 ; J: 3.6/i, and 4.5/1 data points arc 0.02, 0.01, 0.01, 0.01, 0.01, and 0.01 mag, respectively. 



ies bv iDrorv et alJ l)2004b|) . who establish the reliability 
of the "photometric masses" . 

In this paper we extend, for the first time, the com- 
parison between dynamical and photometric masses to 
z = 1. We us e our sample of early-ty pes with veloc- 
ity dispersions ijvan der Wei et aLll2005^ . and rest-frame 
UV-through-IR photometry. We investigate the scat- 
ter between dynamical and photometric masses, both at 
z = and z = 1, and how this depends on the pho- 
tometry included in the SED fits, and the model as- 
sumptions. Furthermore, we investigate systematic dif- 
ferences between dynamical and photometric masses. We 
note that photometric masses are never absolute, as dark 
matter fractions and numbers low-mass s tars are uncon- 
strained from the photometr y (see, e.g.. iGerhard et abl 
I200H iCappellari et alJ 12005(1 . However, photometric 
masses have been used to obtain "relative masses" at 
higher redshif ts, and to determine t he relative mass evo- 
lution (e.g.. iDickinson et alJl2003t IRudnick et al] 12003 
IBell et al.ll2004t IRudnick et al.ll2006ft . We test here ex- 
plicitly whether the relative photometric masses at z = 
and z = 1 are consistent with the dynamical masses. 
The only assumption we make is that high- and low- 
rcdshift early-types have the same kind of stellar popu- 
lation (IMF, metallicity), with (obviously) different ages. 
This assumption underlies all work with relative photo- 
metric masses. 

We specifically address the question whether the rest- 
frame near-infrared (near-IR) helps to constrain the 
masses. Because the near-IR is less sensitive to extinction 
than the optical, extending SED fitting to the near-IR 
helps to lift the degeneracy between age and extinction. 
On the other hand, the stellar population models are less 
reliable at wavel engths longer than lfj, than in the optical 
l)Maraston 2005'). Furthermore, optical-to-near-IR colors 
have been shown to correlate less well with M/L than op- 



tical colors (IBell fc de Jongf2 001). The advance of IRAC 
l|Fazio et al.ll2004|) provides access to the rest-frame IR 
at high redshift, and it is assumed that this will allow 
for 'cheap' mass determinations through SED fitting for 
large numbers of distant galaxies. It is clear that ex- 
tending the fit SEDs to the IR and thus obtaining mass 
estimates needs to be tested. 

This paper is organized as follows: we describe the dy- 
namical masses of the galaxy samples in Section [21 The 
photometry, the stellar population models, and our fit- 
ting method to obtain stellar masses are described in Sec- 
tion[3] In Section^lwe present our results, discussing the 
consistency of the models with our empirical results, and 
the dependence on the fit wavelength range. In Section 
03 we discuss the biases that are revealed by this work, 
and how this affects estimates of high-z galaxy masses 
and the evolution of the mass density. Through out we 
use th e Vega magnitude system (based on the iKuruczl 
(1992) AOV model spectrum), and the concordance cos- 
mology, (n M ,n A ,h) = (0.3,0.7,0.7). 

2. DYNAMICAL MASSES OF EARLY- TYPE GALAXIES 

Dynamical masses are computed as Md yn = 
CreffCc/G, where C is a constant, a c is the central ve- 
locity dispersion, r cr f is th e effective radiu s and 67 is 
the gravitational constant. iKochanekl (fTil!) 1 i . assuming 
spherical symmetry, no rotation and fixed anisotropy, has 
shown that the velocity dispersion of the dark matter in 
elliptical galaxies equals the central line-of-sight veloc- 
ity dispersion, as measured within a 2" x 4" sized aper- 
ture. This implies C = 4 if Md yn is assumed to be twice 
the mass within the effective radius. In terms of the 
dispersions we use (corrected to an 3'.'4-diameter aper- 
ture at the distance of the Coma cluster) , we would have 
to use C = 4.11, taking the differences in aperture and 
distance between the samples into account. In the lit- 




Fig. 1. — Various model spectra for a range of ages and parameters. The labels 'BC03' identify Bruzual & Chariot (2003) 
models; the label 'M05' identifies the Maraston (2005) model. Metallicity and IMF slope are also indicated [x = 1.35 corresponds 
to the Salpeter IMF). The indicated ages run from top to bottom. M05 does not provide realistic spectra for A > 2.5/i. Note 
the difference in optical-to-near-IR slope between the BC03 and M05 models. 



erature, C = 5 is used more often (e.g.. Uoreensen et al.1 
11996(1 . Therefore, we choose to use C — 5 throughout this 
paper in order for the masses to remain consistent with 
previous studies. This number is based on an isotropic 
density model, projected onto an r 1 ' 4 -law and integrated 
out to 2.5rf;ff, which 75% of the total mass, (see, e.g., 
iCappellari et aljl2005(l . In units of solar masses, the dy- 
namical mass becomes M/Mq — 1.17 x 10 6 r c ga^, with 
r e g in kpc and a c in km s _1 . We note, however, that the 
precise value of C does not affect our analysis as long as 
it does not evolve w i th red shift. 

Ivan der Wei et all l)2005fl provide internal velocity dis- 
persions and structural parameters for 29 galaxies in 
the Chandra Deep Field-South (CDF-S). We exclude 
four galaxies with late-type morphologies, four early- 
types with spectra with low signal-to-noise ratios (S/N < 
12 A -1 ), and one early- type that falls outside the IRAC 
mosaic. We then have a sample of 20 early-type galax- 
ies with a median redshift of z — 0.98. The median 
mass of this sample is 1.5 x 10 11 M Q . As a low-redshift 
comparison sample we use 23 early-type galaxies in the 
near by (z = 0.024) Com a cluster with measured a c and 
r eff ^Dressier et a1.l1 987HFa.her et a.l .11 989ft. The median 
mass of this sample, M — 1.8 x lO n M0, is comparable 
to that of the z ~ 1 sample. 

3. DERIVATION OF STELLAR MASSES FROM 
PHOTOMETRY 

3.1. Photometric Data 

A large range of photometric data is available for 
the CDF-S. GOODS provides ACS imaging in 4 filters 
(F435W, F606W, F775W F850 LP, hereafter 6 435 , «606, 
*775j £850 <|Giavalisco et al.l l2004)). ESO provides J- and 
K-haxid imaging, and IRAC GTO observations from the 
Spitzer Space Telescope are available (channels 1-4; 3.6/i, 
4.5/i, 5.8/i, a nd 8.0/x, respective l y). T he photometry is 
described by [van der Wei et alJ <)2006j) and is presented 
here in Tabled The magnitudes are measured within a 
fixed aperture with a diameter of 5'.'0, using registered 
and PSF-matched images. 



For the local sample of Coma galaxies we need photom- 
etry that samples a similar rest-fra me wavelength range 
as we have for the z ~ 1 sample. iFaber et alJ l|1989D 
and iScodeggio et ftT| l|1998 [) provid e the optical surface 
photometry. and lPahre et aiTl|1998fl provide the if -band 
surface photometry. We use the effective surface bright- 
nesses and effective radii in the B-, V-, I-, and Zf-bands 
to compute the colors within apertures of 67" , used by 
Faber et al., which typically corresponds to 4r c g at the 
distance of the Coma cluster. The average colors are 
B -V = 0.96, V - I = 1.18, and I - K = 1.96. The 
Coma galaxies are the only sample in the literature with 
K-band photometry suitable for this study. The fact that 
these are cluster galaxies, and not field galaxies as the 
galaxies in our distant sample are, does not limit the 
interpretation of our results. The differences between 
the two samples are reproduced through modeling their 
SEDs, adopting the relevant parameters which might dif- 
fer (such as age, SFH, and dust content) as free param- 
eters in the fit. The only assumption is that they have 
similar stellar populations in terms of IMF and metallic- 
ity. 

3.2. Stellar Population Models 

For our purpose we need stellar population models that 
provide synthetic spectra ove r a large wavelength range 
(0.3 — 4/i). In recen t years. | Bruzual & Charlotl ( 2003) 
(hereafter, BC03) and lMarastonl l|2005H (hereafter. M05), 
have provided such spectra. M05 provide realistic spec- 
tra only up to 2.5/x, limiting our SED fitting range to 
the 4.5/i-channel. This is not a severe limitation as the 
photometric errors in the two longest-wavelength chan- 
nels are much larger than in the other two channels (see 
Table QJ, and observed 4.5/i at z ~ 1 corresponds to the 
reddest filter (K) available for the nearby sample. 

Significant differences between the models occur 
for all ages of interest for this study. For 
ages 0.5 — 2 Gyr, the M05 model is much red- 
der in the optical-to-nearlR than BC03, because of 
the very different implementation of the Thermally- 
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Pulsating Asymptotic Giant Branch (TP-AGB) phase. 
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Fig. 2. — Evolution of B — I (a) and B — K (b) with age for 
different models and different model parameters. The differ- 
ences between the BC03 and M05 models for similar sets of 
parameters are striking. For ages < 3 Gyr differences up to 
0.5 mag occur in B — I, in B — K even up to f .4 mag. For 
older ages, the agreement is good in B — I, whereas in B — K 
the difference is ~ 0.5 mag. The tracks are limited to ages 
from 0.5-20 Gyr. 

For ages older than 3 Gyr, the M05 model is bluer in 
optical-to-near-IR colors than the BC03 model because of 
the coo ler Red Giant Bran ch (RGB) stars of the Padova 
tracks (|Girardi et alJl200Cl u sed by BC03), with respect 
to those of the Cassisi tracks ijCassisi et all2000l used by 
M05). M05 adopt the fuel consumption approach, and 
calculates luminosity contributions from different stellar 
types by the amount of fuel used during a certain evo- 
lutionary stage, whereas BC03 follow the isochrone syn- 
thesis approach, and characterize the properties of the 
stellar population per mass bin. 

We consider SSP models from M05 with solar and 
super-solar metallicity (Zq and 2.2Z Q ) and a Salpeter 
IMF with mass limits 0.1M Q and lOOAf© (these limits 
are used throughout the paper) . For the BC03 model we 
also consider models with two different metallicities (Zq 
and 2.5Zq) and, in addition, a solar-metallicity model 
with a 'top-heavy', or 'flat' IMF with slope of x — 0.35 
(instead of x = 1.35, which is the Salpeter IMF). The 
model with a flat IMF has been shown to provide a better 
match to the evolution of M/L in the optical an d in the 
near-IR simultaneously <)van der Wei et al.l200rl . We do 
not consider IMFs with different shapes at the low-mass 
end. The reason is that varying this shape changes all 
mass estimates by a constant factor. In other words, the 
high- and low-redshift samples are affected similarly by 
the choice of parametrization of the low-mass end of the 
IMF. 

In Figure^we show the spectra of the three BC03 mod- 



els and the M05 model for ages of 0.5, 1, 2, and 5 Gyr. 
The difference between the BC03 models with different 
metallicities, shown in the left-hand panel, is subtle, ex- 
cept for the longest wavelengths. The BC03 models with 
different IMF slopes (middle panel) show much larger 
variation. This is primarily a difference in overall en- 
ergy output, which is simply due to the larger numbers 
of giants in the case of a flat IMF, especially at young 
ages. A secondary, but actually more interesting differ- 
ence is the different evolution of the slope, i.e., color, of 
the SED in models with different IMFs. Such differences 
in color evolution are also apparent when comparing the 
BC03 models to those of M05, shown in the right-hand 
panel. For ages younger than 2.5 Gyr, the M05 model 
predicts much higher luminosities in the iC-band than the 
BC03 model, whereas the optical luminosities are simi- 
lar. The differences in color evolution are demonstrated 
more clearly in Figure For ages older than 3 Gyr the 
optical colors from different models agree well. For ages 
younger than that differences of up to 1 magnitude in 
B — I occur between BC03 and M05 models with the 
same model parameters. For the near-IR differences oc- 
cur for all ages and can increase up to 1.4 magnitudes in 
B-K. 

3.3. Fitting Method 

We use redshifted model spectra to compute appar- 
ent magnitudes, m mo< j, allowing the age (and, optionally, 
the star-formation history (SFH) and the dust content) 
to vary, and normalizing the calculated magnitudes to 
match the observed magnitudes, m b s , to obtain the pho- 
tometric mass. m i, s are the color magnitudes given in 
Tabled including a correction from fixed aperture mag- 
nitudes to total magnitudes as measured in the if-band. 
The resulting total fluxes are multiplied by 0.75 to match 
the aperture in which the dynamical mass is calculated 
(see Section EJ. The derived stellar mass includes the 
dark, compact remnants of massive stars, but not the 
gas lost due to stellar winds and supernovae ejecta. The 
best-fitting model is selected on the basis of the root 
mean square of m mo d ~ mobs (RMS) , weighing with the 
inverse square of the photometric errors. A certain mini- 
mum error is assumed to avoid that the data points with 
the largest errors are effectively ignored in the SED fits. 
When the IRAC data are used in the fit, a minimum er- 
ror of 0.10 mag is adopted. When the IRAC data are 
not used, we use 0.03 mag. The choices are based on 
the relative uncertainties in the photometric zeropoints 
of the different datasets. However, the precise value of 
the minimum error does not profoundly affect the fitting 
results for the samples as a whole 

For the SSP models described in the previous section, 
there is a unique relation between color and age. Sub- 
sequently, age determines M/L, and thus M p hot- The 
differences between the models, discussed in the previous 
section and shown in Figure |2 will cause mass estimates 
derived from SED to differ. Most notably, the different 
trends with age will cause systematic biases in mass es- 
timates of young (high- 2: ) galaxies relative to old (low-z) 
galaxies. 

In Figure we show, as an example, the fitting re- 
sults for object CDFS-2, using the SSP models. In the 
top panel we show how log(M p ^ ot ) of the best-fitting 
BC03 model (with solar metallicity and a Salpeter IMF) 



van der Wei et al. 
changes if different wavelength ranges are used in the fit. 
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log(M/M Q ) 

Fig. 3. — The quality of the SED fit, quantified by the root 
mean square (RMS) of the difference between the calculated 
model magnitudes and the observed magnitudes, of object 
CDFS-2 at z = 0.96, as an illustration of the fitting method. 
The upper panel shows how the results vary by changing the 
wavelength range included in the fit, using the BC03 model 
with solar metallicity and a Salpeter IMF. dashed line: result 
for fitting 1775, 2850, and J. solid line: result for fitting 1775, 
0850, J, K, 3.6/i and 4.5/i. dotted line: result for fitting 1775, 
£850, J, K, 3.6/i, 4.5/i, 5.8/x and 8.0/i. Including the rest- 
frame near-IR in the fit leads to a somewhat higher M/L. The 
bottom panel shows how the results vary when modeling 1775 
through 4.5/x using different models, solid red line: result for 
the BC03 model with Solar metallicity and a Salpeter IMF. 
dashed red line: result for the BC03 model with super-solar 
metallicity. dotted red line: result for the BC03 model with a 
top-heavy IMF. solid blue line: result for the M05 model with 
solar metallicity and a Salpeter IMF. Different models clearly 
yield different M/L ratios. The broad minimum obtained for 
the M05 model can be explained by the slow evolution of the 
optical-to-near-IR color (see Figure and the consequently 
poorly constrained age. 

In the bottom panel we show how the results vary from 
model to model, fitting the SEDs from 1775 to 4.5/i. 
\og(M p hot) of the best-fitting models vary by a factor 
of 3, indicating the level of the systematic uncertainty 
in the photometric mass estimate. Note that the quality 
of the fits is generally good: the RMS of the best-fitting 
model is typically only ~ 0.02 mag. 

To further illustrate our SED fitting method, we 
show in Figure 01 the SEDs of three z ~ 1 early- 
types and two different, best fitting model spectra for 
each of those. Those model spectra are the SSP 
models with solar metallicity and a Salpeter IMF of 
BC03 and M05, which are fit to the 1775, zgso, and 
J data-points, i.e., the rest-frame optical. While sim- 
ilar in the rest-frame optical, the two models dif- 
fer from each other in the rest-frame near-infrared. 
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Fig. 4. — SEDs and model spectra of three typical z ~ 1 
early-type galaxies. The model spectra, which are for SSPs 
with solar metallicity and a Salpeter IMF, are fit to the three 
shortest wavelength data-points shown here (1775, zsso, and 
J), with age as the only free parameter besides the normaliza- 
tion. The BC03 model under-predicts the fluxes in the redder 
bands, whereas the M05 model makes better predictions. 

The BC03 spectrum predicts systematically lower flux 
levels than the M05 model. The K, 3.6/im, and 4.5/xm 
data-points agree better with the M05 model than with 
the BC03 model. As a consequence, when we fit the 1775 
through 4.5/xm SEDs of these galaxies, the best-fitting 
age increases for the BC03 model in order to match the 
redness of the SEDs. No older ages are found when the 
M05 is used. This behavior is typical for the galaxies in 
our sample. As we will show in the subsequent sections, 
whether or nor including the near-infrared in the SED fits 
has profound consequences on the resulting photometric 
mass estimates. 

4. SYSTEMATIC UNCERTAINTIES IN PHOTOMETRIC 
MASS ESTIMATES 

4.1. Rest-Frame Optical SED Fits 

In this section we compare photometric and dynami- 
cal mass estimates. We begin by fitting the rest-frame 
optical SEDs, as previous work has shown that a single 
rest-frame opti cal color like B — R c an be used well to 
estimate M/L ( Be ll fc de Jondl200lD For the local sam- 
ple we use the B-, V-, and /-bands in the fit, and for the 
distant sample we use 1775, zsso, and J- 

We start with mass estimates resulting from fits of 
the SSP model from BC03 with solar metallicity and 
a Salpeter IMF. The left-hand panel of Figure 03 shows 
Mdyn vs. Mphot for the sample of local early-type galax- 
ies. On the right-hand side we show the same relation 
for the distant galaxies. 

To quantify the difference between Mphot and M^yn, we 
introduce the parameter Si = log(M p f lot / Mdy n ) , where 



6 



Comparing Dynamical and Photometric Masses at z = 1 



12 - 



o 
si 
tx 



2 

\ 

SO 
o 



11 



10 



1 1 1 I 1 1 1 

Coma galaxies 


' ' ' 1 ' 


1 , , , 
Distant 


, 1 , , , , 1 , 
field galaxies // 


~ <z> = 0. 
Ont i or 1 


02 


/ • — 

• 

• ..*' 


~ <z> = 


0.93 //- 
// - 

,/♦ - 




• 




- + 


_ hi?-/** 


•r' 


••'f 
• 




* 


*✓ 


- (a) /" 






-(b) // 


A = -0.11±0.07. 


/ , , , 


(7 = 

, 1 , 


0.17±0.02 _ 
i i i 1 i 


- // 
■V , , , 


<7 = 0.25±0.03 _ 
i 1 i i i i 1 i 


10 


11 


12 


10 


11 12 



log(M/M ) (dyn.) 



Fig. 5. — (a) Comparison between dynamical and photometric masses of early-type galaxies in the Coma cluster at z — 0.024. 
The photometric masses are obtained by fitting the photometric SEDs (B, V, I) by model spectra from BC03 for a SSP with 
solar metallicity and a Salpeter IMF. Age and mass are the only free parameters. The dotted line indicates the best-fitting line 
with slope 1, indicating the average deviation, a is the standard deviation from that line, (b) Comparison between dynamical 
and photometric masses of the distant field early-type galaxies. Photometric masses are obtained by fitting the photometric 
SEDs from 1775, ^sso, and J by the same BC03 model as used for the local galaxies. The dotted line, the same as in panel (a), 
is shown as a reference. The dashed line is the best-fitting line with slope 1 to the distant field sample, excluding two galaxies 
with velocity dispersions lower than 100 km s _1 . The errors on the photometric masses are obtained through Monte-Carlo 
simulations. 



the subscript i = or i = 1 refers to the redshift of the 
sample. Sq = —0.10 ±0.04, which means that photomet- 
ric masses are 20 ±8% smaller than dynamical masses for 
the low-redshift sample. For the high-z sample we find 
5\ = -0.20±0.07. The fact that S t < may be explained 
by dark matter, or an underestimate of the number of 
low-mass stars in these galaxies. It is more relevant that 
Si and $o are consistent with each other: the photome- 
try does not allow us to constrain the low-mass end of 
the IMF and the dark matter fraction. Therefore, we 
normalize the iow-z photometric masses such that, on 
average, they are equal to the dynamical masses. The 
z ~ 1 photometric mass estimates are changed by the 
same amount such that the normalized offset between 
Mphot and Mdyn at 2 ~ 1 can be expressed as A = S\— Sq. 
Thus, for the above, we have A = —0.11 ± 0.07, i.e., 
the normalized z ~ 1 photometric masses are consistent 
with the dynamical masses (they are marginally smaller, 
by 22 ± 14%). The fact that A is consistent with zero 
means that the differences between the SEDs of the low- 
and high-redshift samples are correctly transformed into 
a difference in M/L by this BC03 model. As can be 
seen, M p hot and Md yn correlate well; the scatter around 
the average M p hot/Mdyn is a factor of 1.50 for the local 
sample and a factor of 1.78 for the distant sample (these 
and all other modeling results for the high-z sample are 
summarized in Table 

In order to investigate systematic effects, we perform 
SED fits with different model parameters. We show the 
results for the high-redshift galaxies in Figure |SJ If we 
fit a super-solar metallicity BC03 model (panel a) we 
obtain 81 = -0.32 ± 0.05 and A = 0.00 ± 0.06. If we 



adopt a flat IMF (panel b), we find Si = -0.70 ± 0.06 
and A = -0.37 ± 0.07. Thus, if IMF and metallicity 
are unconstrained the average photometric mass of the 
high-z sample is uncertain by at least a factor of 5 in 
an absolute sense. The differential M/L from z = 1 to 
z = is uncertain by a factor of 2.3. 

Now we explore the M05 models, using the same ages, 
metallicities and IMF as for the BC03 models. In panel c 
of FigurelHlwe show the z ~ 1 results if we adopt the M05 
model with solar metallicity and a Salpeter IMF (see also 
TableEJ). We find that A = -0.23 ±0.06, which is signif- 
icantly smaller than what was found for the BC03 model 
and also significantly smaller than zero: the normalized 
high-z photometric masses are a factor of 1.7±0.2 smaller 
than the dynamical masses. We find better agreement for 
the super-solar metallicity M05 model (panel d of Figure 
EJ): A = -0.08 ±0.08. The disagreement between the 
BC03 and M05 models implies a systematic uncertainty 
of a factor of 1.3 in the normalized high-z mass estimates. 
This is much smaller than the uncertainty due to uncon- 
strained metallicity and IMF. The scatter in M p hot/Mdyn 
is similar to that obtained with the BC03 models (a fac- 
tor of 1.65-1.85). 

To illustrate the results described above, we show in 
Figure [7^ the relation between optical color (B — I) and 
M/L for the models and the observed galaxies. The 
BC03 models with a Salpeter IMF reproduce the col- 
ors and M/L rather well, whereas the BC03 model with 
the flat IMF and the M05 models do not. For the high- 
z galaxies, the M05 model under-predicts M/L, which 
explains that this model gives A < 0. 
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TABLE 2 

SED FITTING RESULTS 



Model 


A 




<5i < rms > 


< Age > 


< SFH > 


< A v > 










mag 


Gyr 


Gyr or % 


mag 






2 = 1: 


«775 , 2g50 , J 










BC03, Zq 


-0.11 ±0.07 


0.25 ± 0.03 


-0.20 ±0.07 


0.03 


2.0 






BC03, Zq, x = 0.35 


-0.37 ±0.07 


0.27 ±0.03 


-0.70 ± 0.06 


0.03 


1.8 






BC03, 2.5Z Q 


0.00 ± 0.06 


0.24 ± 0.02 


-0.32 ±0.05 


0.04 


1.1 






M05, Z e 


-0.23 ±0.06 


0.22 ±0.03 


-0.33 ±0.05 


0.05 


1.7 






M05, 2.2Z e 


-0.08 ±0.08 


0.27 ±0.03 


-0.48 ±0.06 


0.05 


0.7 






2 = 1: i 775 , 2 850 , J , K , 3.6^1 , 4.5/U 


BC03, 


0.15 ±0.06 


0.22 ± 0.04 


0.08 ± 0.05 


0.10 


3.9 






BC03, Zq, t 


0.03 ±0.06 


0.23 ± 0.03 


-0.04 ±0.03 


0.04 


3.0 


0.97 


0.56 


BC03, Z 0) double burst (10%) 


0.07 ± 0.06 


0.23 ± 0.04 


0.00 ± 0.05 


0.07 


2.9 


10% 


3.52* 


BC03, Zq, double burst (30%) 


0.01 ±0.07 


0.24 ±0.03 


-0.06 ±0.04 


0.06 


2.5 


30% 


1.87* 


BC03, Zq, x = 0.35 


-0.02 ± 0.07 


0.22 ± 0.03 


-0.39 ±0.05 


0.09 


3.4 






BC03, 2.5Z Q 


0.23 ± 0.07 


0.25 ± 0.04 


-0.23 ±0.06 


0.08 


1.3 






M05, Zq 


-0.26 ±0.07 


0.24 ± 0.04 


-0.19 ±0.06 


0.09 


3.1 






M05, Zq, r 


-0.24 ±0.07 


0.27 ±0.04 


-0.17 ±0.06 


0.05 


3.7 


0.15 


0.36 


M05, Zq, double burst (10%) 


-0.28 ±0.07 


0.29 ± 0.03 


-0.21 ±0.07 


0.07 


2.2 


6.1% 


1.74* 


M05, Zq, double burst (30%) 


-0.32 ±0.07 


0.27 ±0.03 


-0.25 ±0.06 


0.05 


2.3 


23% 


1.41* 


M05, 2.2Z 


-0.08 ±0.08 


0.27 ±0.05 


-0.37 ±0.06 


0.09 


1.5 










2 = 1: 6435 , 


^606 , *775 , 2850 , J 










BC03, Zq, x = 1.35 


-0.18 ±0.07 


0.28 ± 0.03 


-0.28 ±0.07 


0.30 


1.8 






BC03, t, A v 


-0.05 ± 0.07 


0.24 ±0.03 


-0.15 ±0.06 


0.20 


2.1 


0.29 


0.28 


BC03, double burst, A v 


-0.03 ± 0.07 


0.28 ± 0.03 


-0.12 ±0.06 


0.24 


2.7 


15% 


0.60* 


M05, Zq, x = 1.35 


-0.31 ±0.07 


0.28 ± 0.04 


-0.41 ±0.06 


0.24 


1.4 






M05, r, A v 


-0.23 ±0.07 


0.25 ± 0.03 


-0.33 ±0.06 


0.22 


1.7 


0.13 


0.22 



NOTE. — SED fitting results for the high-z sample. The photometry used in the fits is indicated for the z ~ 1 sample. For the z — sample 
similar rest-frame wavelength coverages arc used (sec text for the exact filter sets). Si is the average log(Af p hot fMd yn ), and eri is the scatter therein. 
A = Si — 5q. Wc use 5q from single burst models to calculate A, also in case of more complicated SFHs and dust-extinction, x — 0.35 indicates 
models with a fiat IMF (with slope 0.35 instead of 1.35 which is the slope of the Salpctcr IMF). SFH indicates cither the average time-scale r of 
the exponentially declining star-formation rate, or the average percentage of stellar mass involved in a secondary burst. This percentage ranges 
between 0% and the adopted burst strength, which indicates that not all galaxies arc necessarily fit better by a double burst model than by a single 
burst model. Ay is the extinction in the rest-frame V band. Those values marked with a * indicate the attenuation of the secondary burst, the 
primary burst is assumed to be dust-free. 
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Fig. 6. — Comparison between dynamical and photometric masses of distant field early-type galaxies as in Figures^] including 
«775, Z850, and J in the fits. The photometric masses are obtained by using different models. Panel (a) shows the results for 
the BC03 model with super-solar metallicity (2.5Z ) and a Salpeter IMF. Panel (b) shows the results for the BC03 model with 
solar metallicity and a top-heavy IMF (with x — 0.35 instead of x = 1.35 (Salpeter)). Panel (c) shows the results from the M05 
model with solar metallicity and a Salpeter IMF. Panel (d) shows the results from the M05 model with 2.2Zq and a Salpeter 
IMF. A is inconsistent with zero in panels a and c, which indicates a discrepancy between the models and the observed SEDs 
and M/L. Also, we note that 8i is smaller than zero in all cases (the dashed and dotted lines lie below the line M p h ot = Md yn ), 
which indicates the presence of dark matter or a larger number of low-mass stars than assumed in the Salpeter IMF. 
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Fig. 7. — The relation between color and M/L for the mod- 
els used to infer photometric mass estimates, compared with 
observed galaxy colors and dynamically obtained M/L. Open 
circles represent the local early-type sample; filled circles rep- 
resent the z ~ 1 early-type sample. Panel (a) shows M/Lb 
and B - 1. The BC03 models with a Salpeter IMF agree bet- 
ter with observations than the M05 models. Panel (b) shows 
M/Lk and B — K. In this case, the M05 model agrees better 
with the data than the BC03 models. 



4.2. Rest-Frame Optical and Near-IR SED Fits 

We now extend the SED fits to the rest-frame K- 
band. The local galaxies are fit with the B-, V-, I- 
, and A"-band photometry. The distant galaxies are 
now fit with ^755 through 4.5/i. Similar to Figure 03 
we show in Figure |H1 the results if we adopt the BC03 
model with solar metallicity and a Salpeter IMF. We 
now find Si = 0.08 ± 0.05, i.e., the high-z photometric 
masses have increased by almost a factor of 2 with re- 
spect to the photometric masses obtained from the opti- 
cal SEDs alone. This is not the case for the local sample 
(So = — 0.07± 0.03), which implies that A has increased 
significantly, to A = 0.15 ± 0.06. Thus, the normalized 
high-z stellar masses are now a factor of 1.4 larger than 
the dynamical masses. A similar increase in A is found 
for the other BC03 SSP models (see panels a and b of 
Figure Hand Table EJ) For the model with a flat IMF 
this means that A is consistent with zero, but we have 
shown in the previous section that such a model provides 
a very poor fit to the rest-frame optical colors and M/L. 

For the M05 SSP models (panels c and d of Figure |§J) 
the effect of including the near-IR is very different. It 
is still true that the high-z photometric masses increase 
somewhat, but that is also the case for the local sample. 
The net effect is that A does not change significantly. For 
the solar metallicity model we find A = —0.26 ±0.07 and 
for the super-solar metallicity we find A = —0.08 ± 0.08. 

The most important consequence of the above results 
is the large difference between the masses inferred from 



the BC03 and M05 models with the same parameters. 
Normalized masses at z ~ 1 obtained with a BC03 model 
are 2-2.5 times larger than normalized masses obtained 
with a M05 model. This implies that there is an intrinsic 
systematic uncertainty in photometric mass estimates of 
at least a factor of 2.5. We stress that our results are 
not caused by the different redshifts of the samples: the 
same rest-frame wavelength is sampled for the low- and 
high-redshift samples. The age difference between the 
samples reveals the systematic problems. In principle, 
tests such as described here do not require a range in 
redshift, but in the local universe there is no suitable 
sample of relatively simple stellar systems with ages of 
1-2 Gyr and dynamically measured masses. 

In Figure 03 we show the relation between B — K and 
M/Lk for the models and the observed galaxies. The 
M05 models reproduce the galaxy colors and M/L better 
than the BC03 models. Note that in the optical the 
reverse is the case (see Figure^)- Qualitatively, we may 
understand that extending SED fits to the near-IR does 
not constrain mass estimates much further: the optical- 
to-near-IR color (e.g., B — K) does not constrain M/L 
very well. Reversely, the M05 model, which predicts the 
lack of strong evolution in B — K cannot distinguish well 
between galaxies with low and high M/L. As a side- 
effect, the scatter in M p hot/Mdy n is similar to that found 
when the SED fits are restricted to the optical, a factor 
of 1.65-1.85 for all models described above. 

The relatively high masses inferred with the BC03 
model for the z ~ 1 galaxies are the result of the rel- 
atively high near-IR luminosities of the galaxies. Ac- 
cording to the model these imply old age and high 
M/L (see Figures |2 and EJ) . Another possible explana- 
tion for the high near-IR luminosities is the presence of 
dust-obscured stellar populations that do not contribute 
to the optical luminosity. Generally speaking, models 
with more complex star-formation histories might yield 
younger ages and lower M/L. We explore models with 
exponentially declining star-formation rates and models 
with a dust-obscured, secondary burst of star-formation 
that occurs 2 Gyr after the initial burst. 

We adopt a fixed contribution to the stellar mass of the 
secondary burst and leave age and extinction as free pa- 
rameters. Leaving the secondary burst strength as a free 
parameter results in a very large scatter in M p hot/Mdyn 
as the three free parameters are degenerate. If we assume 
that a secondary burst accounts for 10% of the total stel- 
lar mass, we find A = 0.07 ±0.06 (where Sq is taken from 
the results obtained with the SSP model as the effect of a 
secondary burst does not affect the SEDs of the old local 
early-types) . The galaxies responsible for the decrease in 
A typically have a 0.1 Gyr old secondary stellar popu- 
lation that is highly obscured (Ay ~ 4). Similar results 
are obtained if the secondary burst strength is increased 
to 30% (A = 0.01 ± 0.07). Models with exponentially 
declining star-formation rates also improve the results 
for the BC03 model if we allow for the presence of dust 
(A = 0.03 ± 0.06). Hence, models with a more com- 
plex SFH can successfully reproduce the SEDs of z ~ 1 
early-type galaxies that yield photometric masses that 
are consistent with their dynamical masses. Also, the 
quality of the fits is better for the models with more com- 
plicated star-formation histories than for the SSP model 
(see Table EJl. We note that these type of models do not 
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Fig. 8. — A comparison between dynamical and photometric masses as in Figure^ but now including the rest-frame near-IR 
in the SED fits. For the local Coma sample, the B-, V-, I-, and K-band data-points are used. The average photometric mass 
is not affected by the changed wavelength coverage: the dotted line has not shifted significantly with respect to Figure For 
the distant sample, 4775 through 4.5/i are used to obtain photometric masses. There is a systematic increase with respect to the 
masses obtained with the optical SEDs alone: the dashed line has shifted with respect to Figure Now, A > 0, which implies 
that stellar masses at z ~ 1 are overestimated relative to z — 0. 



improve the results for the super-solar metallicity BC03 
model (A > 0.2, see Table 0. The M05 model fits do 
not change significantly if a model with a complex star- 
formation history is used (A < —0.2, see Table 

4.3. Rest-Frame UV and Optical SED Fits: Constraints 
on the Star- Formation History 

Finally, we explore the effect of extending the analy- 
sis to the rest-frame UV. The UV is strongly affected 
by (small traces of) star-formation, and is expected to 
cause less accurate mass estimates. On the other hand, 
this sensitivity to young stars may allow to constrain the 
SFH and dust content, which we below adopt as free 
parameters. 

As a baseline, we adopt the BC03 and M05 SSP models 
with solar metallicity and a Salpeter IMF. We now fit the 
SEDs using the 6435, «606j 2775, z 850, and J photometry. 
Because of the problems detailed in the previous section, 
we omit the rest-frame near-IR. The results are shown 
in the left-hand panels of Figure fTOl (BC03) and Figure 
ITTl (M05). The scatter in M p hot/Mdy n is a factor of ~ 
1.9, similar to the scatter obtained without the 6435 and 
t>606 photometry. The quality of the fits has decreased 
significantly: for the BC03 model, the average RMS is 
now 0.30 mag (it is 0.03 mag without the two shortest- 
wavelength filters); for the M05 model this is 0.24 mag, 
instead of 0.05 mag. This indicates that, at least for the 
UV part of the SEDs, an SSP model provides poor fits. 

Now we allow for an exponentially declining star- 
formation rate with timescale r, instead of assuming a 
single burst. Since star- formation is usually related to 
extinction by dust, we consequently adopt Ay as a free 
parameter as well. The results are shown in the middle 
and panel of Figure ITUl fBC03) and the right-hand panel 
of FiguretTTlCM05V The scatter in Mphot/Mdyn decreases 



marginally and is now comparable (a factor of ~ 1.7) to 
what we found when we excluded the UV and assumed a 
dust-free SSP. Also, the fits are somewhat better, as the 
average RMS is now ~ 0.2 mag for both models. The 
fact that the RMS remains much higher than for the op- 
tical SED fits is most likely caused by the large errors in 
the 6435-band photometry. 

In Table El we show the average values of age, r, and 
Ay. The inferred values for the extinction are quite large: 
Ay ~ 0.25 for both models. Age and r are ~ 2.1 and 
~ 0.29 Gyr for the BC03 model with dust and ~ 1.7 and 
~ 0.13 Gyr for the M05 model with dust. This implies 
an average SFR of ~ 1 M Q yr~ x . Assuming that this 
continues to decline exponentially between the epoch of 
observation and the present, the average stellar mass in- 
creases by no more than ~ 0.2% after z ~ 1. The galaxy 
with the highest SFR would increase its stellar mass by 
8.0%. The models imply that residual star-formation is 
not very relevant for the galaxies in our sample, and, 
barring interactions or mergers, have assembled most of 
their stellar mass by z ~ 1. 

Finally, we test a BC03 model in which a dusty, sec- 
ondary burst of star-formation occurs 2 Gyr after the 
dust-free initial burst. This is the same model as used 
in the previous section, but now we adopt the strength 
of the secondary burst as a free parameter as well, along 
with age and Ay: we allow the secondary burst to vary 
in strength between 0% to 30% of the final stellar mass. 
The results are shown in the right-hand panel of Figure 
1101 The scatter in Mphot/Mdyn is a factor of 1.9, which is 
similar to what was found when modeling the rest-frame 
optical SEDs with a dust-free SSP model. 

To investigate the applicability of the double burst 
model, it is more interesting to look at individual galax- 
ies. The fits of 11 out of 20 galaxies clearly improve by 
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Fig. 9. — A comparison between dynamical and photometric masses as in Figure |H1 (fitting rest-frame optical through near-IR 
SEDs) . The same set of models as in Figure |S] are used. For any BC03 model, A increases significantly when the near-IR is 
included in the SED fits. For the M05 model A does not significantly change. 
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Fig. 10. — Comparison between dynamical and photometric masses (using the BC03 model) of distant field early-type galaxies 
as in Figure [3] Here we modeled the rest-frame UV+optical SEDs: all photometric data points between 6435 and J are used in 
the fits. In panel (a) we assume, as a starting point, a dust-free SSP with solar metallicity and a Salpeter IMF. In panel (b), r 
and Ay are free parameters. In panel (c), the strength of a secondary burst that occurs 2 Gyr after the initial burst and Ay 
(only applied to the secondary burst) are free parameters. See Table|5|for the inferred ages, and average values of r and Ay. 



allowing a dusty, secondary burst. With a single burst, 
these seven galaxies are assigned ages of 1.2 Gyr. These 
galaxies are fit significantly better by a ~ 0.3 Gyr old 
secondary burst involving, on average, 12% of the final 
stellar mass and an extinction of Ay ~ 0.65. This young 
population is superimposed on a 2 Gyr older, dust-free 
population. 

5. SUMMARY AND DISCUSSION 

In this paper we have compared photometric and dy- 
namical masses of early-type galaxies at both high [z ~ 
1) and low redshifts. The uncertainties were analyzed 
quantitatively. We start this discussion with summariz- 
ing our main conclusions. 

1) After allowing for a systematic offset, the scatter in 
Mphot/Mdyn, which quantifies the random uncertainty 
in photometric mass estimates of individual early-type 
galaxies at z ~ 1, is a factor of ~ 1.75 (a = 0.25). We 
find no significantly different random uncertainties for 
different models and different wavelength ranges used in 
the fits, i.e., including IRAC data in the SED fits does 
not provide significantly more accurate mass estimates. 



2) For SSP models from BC03 and rest-frame optical 
SED fits we find that the photometric masses of the z ~ 1 
galaxies are consistent with the dynamical masses after 
normalization with respect to the z = galaxies. This 
normalization is the difference between the photometric 
and dynamical masses of the low-z comparison sample. 
When the SED fits are extended to the rest-frame near- 
IR the normalized z ~ 1 mass estimates increase by a 
factor of 2: the z ~ 1 photometric masses are are larger 
than the dynamical masses. For the SSP model with a 
top-heavy IMF (x = 0.35 instead of x = 1.35) there is 
no such discrepancy, but in that case optical SED fits 
underestimate the z ~ 1 masses by a factor of 2.3. 

If we allow for more complex star-formation his- 
tories and dust-extinction, the z ~ 1 mass esti- 
mates decrease, such that the discrepancy seen for 
the BC03 models with a Salpeter IMF and solar 
metallicity disappears. For example, a model with 
a secondary, obscured burst of star formation pro- 
vides good results and increases the quality of the fits. 
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Fig. 11. — Comparison between dynamical and photomet- 
ric masses (using the M05 model) of distant field early-type 
galaxies as in Figure Q>] modeling the rest-frame UV+optical 
SEDs: all photometric data points between 6435 and J are 
used in the fits. In panel (a) we assume, as a starting point, 
a dust-free SSP. In panel (b), r and Ay are free parameters. 
See Table|5|for the best-fitting model properties. 

If correct, these results imply the presence of a signif- 
icant population (at least 10% in mass) of young stars 
(~ 0.1 Gyr) that are highly obscured (Ay ~ 4) in roughly 
half of the z ~ 1 population. It remains to be seen 
whether this is the case. Mid-IR observations will con- 
strain this scenario. 

3) We find very different results when we use the M05 
model with solar metallicity. That model produces nor- 
malized mass estimates from the optical SEDs of the 
z ~ 1 galaxies that are too low by a factor of 1.7. The 
most striking difference with the BC03 model is revealed 
when the SED fits are extended to the near-IR: the nor- 
malized mass estimates at z ~ 1 do not change signif- 
icantly, which implies that the normalized high-z mass 
estimates remain too low. More complex models do not 
change these results significantly. Most importantly, the 
differences between the results obtained with the BC03 
and M05 SSP models with identical model parameters 
imply a systematic uncertainty of a factor of 2.5 in pho- 
tometric mass estimates obtained from near-IR SED fits. 
Hence, extending SED fits from the optical to the near-IR 
does not provide better mass estimates of high- z early- 
type galaxies, because of uncertainties intrinsic to the 
stellar population models. 

4) Adopting the SFH and dust-extinction as free pa- 
rameters, and extending the SED fits to the rest-frame 
UV, do not increase the quality of the SED fits, nor 
decrease the random uncertainty in the photometric 
masses. For our sample, fitting the rest-frame opti- 
cal SEDs is optimal to constrain the masses. Fits of 
UV+optical SEDs of several individual galaxies, how- 
ever, improve significantly by adopting the SFH as a free 
parameter. 

The above results demonstrate that both models have 
problems in either the optical (M05) or the near-IR 



(BC03, barring significant, obscured stellar populations). 
This conclusion relies on two assumptions. First, we 
assume that all local and distant galaxies have simi- 
lar dynamical structures. In other words, we assume 
that invoking 5r c ffcr?/G as Md yn for both low- and high- 
rcdshift early- types does not introduce systematic effects. 
This assumption appears to be reasonable, as FP stud- 
ies have shown that all nearby elli pticals and SOs to fol- 
low the same relat i on (e .g. IJorgensen et al]li.996f) . and 
iTreu &: Koopmansl 1)20041) showed that high redshift el- 
lipticals have the same structure as low redshift ellipti- 
cals. Second, we assume that local and distant galax- 
ies can be fit with models with the same metallicity and 
IMF. By that we ignore, for example, the possibility that 
the galaxies in our z ~ 1 sample evolve strongly in metal- 
licity. However, in order to remove the inconsistency 
found for the BC03 model, the distant galaxies would 
have to be more metal-rich than the local galaxie s . This 
seems unlikely. On the contrary, IJdrgensen et alJ (|2005f) 
show evidence for metal enrichment between z ~ 0.8 and 
the present. We also note that mergers between early- 
type galaxies, which are likely very important in shaping 
the mass- function of th e early- type galaxy population 
(e.g.. Ivan Dokkumll2005|) . do not affect our conclusions 
much because the stellar populations themselves do not 
change. 

The differences between the models are not restricted 
to the range of ages of the galaxies described in this pa- 
per (> 1 Gyr). On the contrary, the largest differences 
are found for even younger ages (i.e., low M/L, see Fig- 
ures and 0). This implies that mass estimates of young 
or star-forming galaxies may also be systematically un- 
certain. This is in particular relevant to studies at z > 2, 
where SED fits are the o nly available method for esti- 
mating galaxy masses (e.g..lForster Schreiber et alJ 2004|: 
I Daddi et all2004lShaplev et al.l2005ULabbe et all2003: 
Moba sher et all 12005^ . This work suggests that more 
than one type of model should be used when perform- 
ing SED fits, as the difference will give an indication of 
the systematic uncertainty. 

It is essential that the models are improved, hopefully 
leading to consistent predictions of the SEDs, given the 
IMF, star formation history and metallicity. Therefore, a 
better understanding of the near-IR properties of stellar 
populations is required. Our work provides a new tool to 
verify the predictions of the models, a method that can 
be exploited more thoroughly by obtaining a large sample 
of early-type galaxies with a wide range in redshift, ac- 
curate dynamical masses and abundance measurements 
from spectroscopy, star-formation activity from UV and 
mid-IR photometry, and mass estimates derived from 
multi- wavelength imaging, including the rest-frame near- 
IR. Only when the model discrepancies are resolved can 
the full potential of IRAC for measuring stellar masses 
of distant galaxies be realized. 
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